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Magnetic resonance techniques provide a powerful tool for controlling spin systems, with applications rang- 
ing from quantum information processing to medical imaging. Nevertheless, the behavior of a spin system 
under strong excitation remains a rich dynamical problem. In this paper, we examine spin resonance of the 
nitrogen- vacancy center in diamond under conditions outside the regime where the usual rotating wave approx- 
imation applies, focusing on effects of multifrequency excitation and excitation with orientation parallel to the 
spin quantization axis. Strong-field phenomena such as multiphoton transitions and coherent destruction of tun- 
neling are observed in the spectra and analyzed via numerical and analytic theory. In addition to illustrating 
the response of a spin system to strong multifrequency excitation, these observations may inform techniques for 
manipulating electron-nuclear spin quantum registers. 

PACS numbers: 



I. INTRODUCTION 

Control over individual electronic and nuclear spins in 
the solid state offers promise for applications in metrol- 
ogy and quantum information science CHS). Manipulation 
of such systems often utilizes a combination of microwave 
(MW) excitation to drive electronic spin transitions and radio- 
frequency (RF) excitation to control nuclear spins O 0171 El • 
Experimental techniques typically rely on sequential RF and 
MW pulses (21 13, but a better understanding of simultane- 
ous MW and RF excitation may enable greater parallelism 
in manipulating few- spin systems. Moreover, multifrequency 
excitation can lead to a variety complex behaviors of funda- 
mental interest |9|. At the same time, control techniques us- 
ing strong-field resonance phenomena have attracted attention 
in systems ranging from trapped atoms ifTOl to superconduct- 
ing devices iTTTTl as well as single spins in the solid state fl2l . 
As transition rates approach excitation frequencies, complex 
dynamics ensue, and theoretical techniques beyond the usual 
rotating wave approximation are required to interpret the sys- 
tem evolution [ 13 ]. In this paper we consider a combination 
of these effects by examining multifrequency excitation of a 
single spin system in a regime where strong-field phenomena 
occur. 

We approach this problem from the point of view of re- 
cent experiments on the nitrogen- vacancy (NV) center in di- 
amond CHS. In particular, we consider a situation directly 
relevant to the study of electron and nuclear spin resonance 
in diamond: simultaneous excitation by weak MW intended 
to drive electronic spin transitions and strong RF intended to 
drive nuclear spin resonance Q. Our experimental observa- 
tions are dominated by the response of the electronic spin, and 
exhibit features characteristic of strong-field excitation such 
as multiphoton resonances and coherent destruction of tun- 
nelling. By comparison to numerical simulations and analytic 
theory, we determine that the behavior we observe involves 
both transverse and longitudinal components of the magnetic 
fields driving the NV spin. Beyond explaining the observed 
resonances, the identification of multiaxis coupling also sug- 
gests techniques for controlling the NV spin system that go 
beyond the usual transverse spin-resonance model. 



II. MULTIFREQUENCY EXCITATION OF THE NV 
CENTER IN DIAMOND 

The electronic spin associated with NV center in diamond 
has recently emerged as a promising candidate for quan- 
tum information processing (T]-[3]| and sensitive magnetome- 
try lfT4HT6l . Its diamond host offers a stiff, nearly spinless 
environment, leading to extremely long electronic spin coher- 
ence times in the range of a millisecond IT71 . Moreover, the 
optical transitions of the NV center allow a high degree of 
spin polarization at room temperature via optical pumping, 
and enable spin-dependent fluorescence for optical spin de- 
tection |[T8ll . With several possible nuclear spins in its envi- 
ronment - 14 N or 15 N and isotopic impurity 13 Cs - the NV 
center offers the opportunity to prepare, manipulate, and mea- 
sure a single NMR molecule in the solid state. 

The ground state of the NV center has an electronic spin 
S — 1 with a natural spin quantization axis S z provided by 
the substitutional nitrogen and adjacent vacancy that comprise 
the defect. In the absence of a magnetic field, the ground state 
m s = spin sublevel lies 2.87 GHz below the m s = ±1 lev- 
els; in the presence of a magnetic field, the NV spin gyromag- 
netic ratio is extremely close to that of the bare electron. The 
optically excited state of the NV defect likewise has an S = 1 
configuration with the same quantization axis and similar gy- 
romagnetic ratio, but its room temperature zero-field splitting 
is only 1.42 GHz l20l (see Fig. la). While electron spin res- 
onance (ESR) has been observed in the excited state l8ll20l, 
the fast optical decay time r ~ 12 ns 1 18 ] limits the resolution 
of such experiments, and the experiments discussed below in- 
volve spin resonance in the ground state configuration. 



A. Experimental observations 

In our experiments, single NV centers are isolated using a 
home-built confocal microscope operating at 532nm for ex- 
citation and 640-800nm for detection. Helmholtz coils or a 
permanent magnet are used to apply static magnetic fields to 
the sample, while oscillatory fields are coupled onto the dia- 
mond by a 25/im copper wire soldered across its surface. We 



2 




MW frequency (GHz) 

FIG. 1: (a) The spin levels of the electronic ground state and opti- 
cally excited state of the NV center. Note that the orbital doublet 
in the excited state is averaged at room temperature where we con- 
duct our experiments 1 19]. The excited state level anticrossing (ES- 
LAC) near 510G is indicated by a dashed oval, (b) The geometry 
of our experiments, illustrating the orientation of oscillating mag- 
netic fields, (c) Multifrequency continuous- wave excitation in low 
magnetic field. Inset shows zoom of boxed region. By extrapo- 
lating from Rabi nutation frequencies at higher MW intensity, we 
estimate the MW electron spin Rabi frequency for this data to be 
Qmw — (27r)0.46 MHz. The RF intensity is significantly higher; 
based on quasistatic analysis (below), we estimate the RF longitudi- 
nal coupling to be Q rf ~ (27r)6 MHz. Each horizontal sweep in 
MW frequency is normalized by its average value, compensating for 
slow fluctuations in fluorescence intensity owing to sample drift as 
the data is collected over several days. 



image an NV center a few fim below the (1 1 1) surface of the 
diamond, choosing an NV center whose axis is perpendicular 
to the surface. Currents through the 25 fim copper wire create 
oscillating MW and RF magnetic fields oriented at an angle to 
the NV axis as illustrated in Fig lb. Furthermore, because of 
the proximity of the narrow wire to the NV center, the ampli- 
tude of the applied oscillatory magnetic fields can exceed 10 
Gauss. 

We examine the NV center spectrum in a multifrequency 
parameter regime relevant for manipulation of coupled elec- 
tronic and nuclear spin systems. Because hyperfine constants 
for nuclear spins in close proximity to the NV center are typi- 
cally in the range of a few MHz (e.g. 2.2 MHz for 14 N, 3 MHz 
for 15 N EH), a combination of few- GHz and few-MHz exci- 
tation is required to resonantly drive both spins (3] [7). Owing 
to tiny nuclear gyromagnetic ratios, the RF oscillatory mag- 
netic field must be quite large (many Gauss) to enable nuclear 
spin transitions on /is timescales, whereas the MW fields driv- 
ing electronic spin transitions may be 10 to 100 times weaker 
to selectively drive individual hyperfine lines. We examine 
simultaneous multifrequency illumination in these frequency 



and power regimes, and find that it produces a complex set of 
resonances dominated by the electron spin response. 

Fig. lc shows a typical multifrequency spectrum for an NV 
center in diamond at low magnetic field, wherein an NV cen- 
ter in a magnetic field of 45 G \ \z is illuminated by 532nm 
light while MW and RF currents oscillate in the wire above it. 
The fluorescence intensity (in the color scale, darker indicates 
lower fluorescence) is recorded as a function of RF and MW 
frequency (urf and uomw)\ darker regions indicate popula- 
tion driven from the m s = state into the m s = — 1 state. 
In the absence of RF excitation, the electron spin resonance 
(ESR) spectrum at this magnetic field would show three hy- 
perfine lines from the three spin projections of the 1 = 1 14 N 
nuclear spin, centered around 2.737 GHz. With RF excita- 
tion, three significant features appear: At low RF frequencies, 
the spectrum is split by 24 MHz; for uj rf > (2tt)5 MHz, 
multiphoton phenomena appear; at moderate RF frequencies, 
the structure becomes complex (see inset) owing to overlap of 
hyperfine lines, multiphoton transitions, and "missing" reso- 
nances such as the one that occurs near ujrf = (27r)5.5 MHz, 
ujmw (2tt) 2.737GHz. 

The multifrequency spectrum of Fig. lc can be simplified 
by polarizing the 14 N nuclear spin. Fig. 2 shows data taken 
at a magnetic field of ~ 510 G \\z where the electron spin 
states show an anticrossing in the excited state (see circled 
region in Fig. la). At the excited state level anticrossing (ES- 
LAC), hyperfine interactions in the excited state of the NV 
center enable resonant electron-nuclear spin flip-flops, effec- 
tively transferring the electron spin polarization to the 14 N 
nuclear spin in a two-step process EH). As a result, in the 
absence of RF excitation, only a single line would be visi- 
ble in the ESR spectrum at 1.45 GHz. While working at the 
ESLAC simplifies the spectrum by eliminating multiple hy- 
perfine lines, it introduces other complications. Because low- 
frequency electron spin transitions are possible in the excited 
state, we avoid applying RF fields during optical illumination; 
instead, MW and RF excitation is pulsed prior to fluorescence 
detection (see pulse sequence in Fig. 2a). The resulting data 
set (Fig. 2a) exhibits the same qualitative features as the low- 
field data, and the triangular pattern of "missing resonances" 
is more pronounced. 

Multifrequency excitation is commonly employed to ob- 
serve double resonance phenomena (22]|23]l, and we note in 
passing that nuclear magnetic resonance (NMR) and electron- 
nuclear double resonance (ENDOR) features can arise in this 
system when ujrf is close to the ground state 14 N nuclear spin 
transition mj = +1 — » raj = (see Fig. 2b). At the ESLAC, 
nuclear spin flips directly correlate to reduced fluorescence, 
creating the horizontal feature observed in Fig. 2b at the NMR 
frequency ujrf ~ (27r)5.1 MHz (71 El- This feature is partic- 
ularly dark because the 10/is pulse is close to a ir pulse on the 
nuclear spin transition, whereas other (time-averaged) transi- 
tions only depopulate the ground state by 50%. Moreover, for 
ujrf near the NMR transition, the raj = state enhances Ra- 
man transitions between the m s = 0, raj = +1 —> m s = 
— 1, raj = states (see Fig. 2b inset). Such transitions occur, 
for example, when ujmw = A—urf — A, where A is the bare 
ESR frequency and A = 2.2 MHz is the 14 N hyperfine split- 
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FIG. 2: (a) Multifrequency spectrum at the ESLAC. Q RF « (27r)15 
MHz and Qmw ~ (27r)0.43 MHz are pulsed for 5 fis prior to fluo- 
rescence detection (sequence shown). Each horizontal sweep in MW 
frequency is normalized by its average value and shifted in frequency 
(using ESR data collected at regular intervals during data acquisition) 
to compensate for drifts in fluorescence intensity and magnetic field. 
Inset shows data superposed with the lines S = tlujrf (red) and ze- 
roes of Js/co RF {^Qrf/urf) (white) for S = ujmw - (27r) 1.4506 
GHz and Q RF = (2tt)15.4 MHz. (b) NMR and ENDOR at the ES- 
LAC. Higher intensities Q M w « (2tt)2 MHz, Q rf w (2tt)20MHz 
and longer pulses (10/is) allow NMR and ENDOR via the mecha- 
nism illustrated in the inset. Examples of weak ENDOR resonances 
are marked with arrows. Note that the ENI 525LA RF amplifier used 
in this experiment does not have a flat response below 2 MHz. 



ting 12H . Indeed, the data exhibits weak ENDOR features in 
the vicinity of ujrf ~ 5 MHz that are roughly 2.2 MHz lower 
in frequency than each multiphoton resonance (see Fig. 2b). 

NMR and ENDOR resonances in the NV center have 
been previously examined at the ground state level-anti- 
crossing 1 22 , 23 ], and they play only a minor role in the over- 
all features observed in the multifrequency data. In particular, 
the three features we consider - the low-frequency splitting, 
multiphoton resonances, and "missing resonances" are not as- 
sociated with nuclear spin transitions. These features of our 
data involve only the electronic spin transitions, and thus re- 
sult from dynamics of an effective two-level system. 



B. Numerical simulations 

To elucidate the origins of the features in the RF/MW spec- 
tra, we numerically simulate the dynamics of a spin subject 
to multifrequency excitation in a parameter regime similar to 
our experiments. Specifically, we consider the evolution of 
the two-level system formed by the m s = and m s = — 1 
spin sublevels governed by the Hamiltonian 

H = ^c^+cr-(f& M w cos (uj M wt) + ^rf cos {uJ RF t + <j>)) 

(1) 

where H = 1, a = {cr x , a y , a z } represents the usual Pauli 
matrices for the pseudo-spin- 1/2 system and ^mw and Orf 
are vectors parallel to the magnetic field orientation for the 
MW and RF excitation respectively. Here, factors of 1/2 are 
chosen so that the transverse orientations of VLmw and Qrf 
correspond to the angular frequency for resonant Rabi nuta- 
tions. 

Experimentally, the MW and RF magnetic fields are ori- 
ented at an angle to the NV axis, but our numerical simu- 
lations indicate that only one component of each vector is 
responsible for the observed resonances. Fig. 3 shows cal- 
culations of the time-averaged population in m s =0 for a 
spin initially polarized into m s = and evolving according 
to Eq. [T]for three different orientations of ^mw and ^rf- 
In particular, the features we observe in our experiments arise 
solely from the x component of ^mw and the z component 

Of f^RF- 

The dependence on MW and RF orientation can be ex- 
plained qualitatively by considering the processes required 
to flip the spin by absorption of one MW and one RF pho- 
ton l25l . Absorption of a photon polarized transverse to 
the NV axis must be accompanied by a spin flip; longitu- 
dinal photon absorption can be thought of as shifting en- 
ergy without changing the spin state. When both MW and 
RF fields are transverse, two-photon processes cannot drive 
a spin transition; only odd-photon-number processes are al- 
lowed, and they are strongly suppressed by the large de- 
tuning of their intermediate states. As a result, Qmw\\x, 
Qrf\\x only significantly drives the single-photon transition 
- a weak three-photon process is also visible in Fig. 3c. 
When allowed by selection rules (i.e. for VLmw \ \x, Qrf | \z or 
&mw \\z, Qrf I \x) the two-photon transition rate will roughly 
scale as VLmw^rf / D where D is the detuning of the in- 
termediate state. For VLmw \ \z, Qrf\ \x, D ~ ^mw, and 
multiphoton transitions are strongly suppressed (see Fig. 3b). 
Only for VLmw \ \&, Qrf\ \z is the detuning D « uorf small 
enough that strong two-photon (and higher) transitions occur 
(see Fig. 3a). Qualitative considerations thus substantiate the 
numerical simulations, indicating that the transverse MW and 
longitudinal RF magnetic field components are responsible 
for the features in our experimental data. 

C. Theory 

While numerical techniques can confirm the behaviors of 
the system and help eludicate their origin, an analytic ex- 
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FIG. 3: Numerical simulation of continous-wave multifrequency 
spectroscopy. All simulations are performed using a fourth-order 
Runge-Kutta integration with a time step of 0. 1 ns over a total time 
of 100/is. EH In each case, ft M w = (2tt)0.5 MHz, Ct RF = (2tt)3 
MHz, and A = (27r)100 MHz, while the simulations average over 
time and <j). The microwave frequency is indicated as a detuning 
from the resonant frequency of 100 MHz; this resonance frequency 
is far lower than the experimental resonance at ~ few GHz, but the 
simulations are not sensitive to this parameter because the RWA is 
valid for the MW excitation. The average population in m s — (for 
a spin initially polarized into m s — 0) is shown in color scale, with 
darker color indicating lower population, (a) SImw \\x and JIrf \\z. 
(b) f&MwII^ and fiRF||£. (c) f&Mw||£and Orf||x. Note that no 
spin transitions occur if both oscillating magnetic fields are oriented 
along z. 



planation provides physical intuition for the features in the 
spectrum and may facilitate development of applications. In 
particular, we analyze the large splitting in the ESR lines at 
low RF frequency, the multiphoton transitions at high RF fre- 
quency, and lastly the "missing resonances" that appear for 
moderate RF frequencies. By use of quasistatic approxima- 
tions and two-level system analysis in the interaction picture, 
we obtain a straightforward explanation of these phenomena. 

The features we seek to understand arise from specific ori- 
entations of the MW and RF field, so we begin by restricting 
Eq.[T]to MW excitation perpendicular to the spin axis (trans- 
verse) and RF excitation along the spin axis (longitudinal). 
Further simplification can be made by applying the RWA to 
the microwave excitation. Choosing a coordinate system to 
eliminate (j>MW and moving to a frame rotating with the mi- 
crowave frequency yields 



H ; 



n 



MW 



a x + tt RF cos (u RF t + (j))(j z , (2) 



where 8 = A — ujmw- The corrections (of order ( Qmw \ ) 
are negligible for the parameter regimes we consider, where 
umw ^> Qrf ~ urf > Qmw- For example, no differ- 
ence is observed for numerical simulations that directly inte- 
grate Eq. [T] and those use Eq. [2] for the parameters given in 



Fig. 3a. Moreover, this Hamiltonian is now of a form that has 
recently received some attention because of its applications to 
superconducting qubits HU [13) : the off-diagonal couplings 
are constant and the diagonal terms are modulated in time. 
Indeed, many of our observations are similar to dynamics of 
microwave driven superconducting devices l26li . 

The behavior at very low RF frequencies in Fig. lc and 
Fig. 2a is readily obtained from Eq. [2] via a quasi-static pic- 
ture of the RF field. In our experiments, the spin population 
evolves coherently during a time t max ~ few [is set by the 
smaller of the dephasing time T 2 * of the electron spin and the 
duration of the excitation pulses. For ujrf <C l/t max , each 
experiment occurs under conditions where the energy split- 
ting 5 + 2Qrf cos {uoRpt + (j)) is approximately constant with 
a random value ranging between S — 2ft rf and S + 2ft rf. 
Owing to the shape of the cosine, the extremal values are 
most probable, and the observed ESR spectrum appears split 
by 2ft rf. Specifically, the probability of a given detuning x 
is proportional to 1 / ^/ 4Tt 2 RF — (x — 5) 2 , which, when con- 
volved with a Lorentzian, is in good qualitative agreement 
with the line shape (see Fig. 4a). Fig. 4b confirms that the 
observed splitting scales with the Rabi frequency Qrf'. solid 
lines illustrate the expected functional dependence Qrf oc 
Vl0 dB / 10 . This quasi-static splitting provides a convenient 
mechanism for characterizing the z component of the RF field, 
and enables the estimates of Qrf we provide for our data. 
Note that these quasi-static effects do not show up in the nu- 
merical simulation because dephasing is not included and the 
timescale for integration far exceeds IIujrf for all calculated 
values of the RF frequency. 

Multiphoton transitions can be readily explained via Flo- 
quet theory l27l . but the parameter regime we investigate in- 
vites an alternate approach (TTJ [13) for analyzing the multi- 
photon and missing resonances. These features occur when 
^mw is the smallest parameter in our system, suggesting a 
natural division of the Hamiltonian into H = H z + Hmw> 
where Hmw is a weak perturbation. The evolution arising 
from the diagonal components H z can be found exactly, al- 
lowing us to move to the interaction picture, where the Hamil- 
tonian becomes 



H 



MW 





a i(nuj RF - 







-i(nujRF—5)t 



(3) 



with A n = ^f^J n (f^ 1 ) (see e.g. Q3) for a full deriva- 
tion). In this picture, the basis states of the interaction 
Hamiltonian are spin eigenstates with time-varying phases 
ex —ifH z dt' 9 transitions induced in the interaction picture 
thus correspond to spin flips in the Schrodinger picture. 

The interaction picture Hamiltonian (Eq. [3]) can be further 
simplified by again applying the rotating wave approxima- 
tion. For weak MW excitation, the off-diagonal components 
of Eq. [3] are much smaller than ujrf, so only the term in the 
sum with S — uujrf « will lead to significant transition 
probability (the others will average to zero over the timescale 
of system response). The different resonances 5 = uuurf thus 
correspond to the multiphoton transitions we observe in our 
data. Furthermore, in the RWA each n-RF photon transition is 
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The spectra seen in our experiments closely match the pre- 
dicted form for CDT resonances. Fig. 4c shows the central 
ESR transition (S = 0) from the data of Fig. 2a, as com- 
pared to the functional form of the expected Rabi frequency 
QmwJo (2Qrf/urf) and numerical integration of the Bloch 
equations for Eq. [2] The maxima in fluorescence (indicat- 
ing suppression of ESR transitions) coincide with the Bessel 
function zeros for Qrf = (27r)15.4 MHz, and the shape of 
the curve fits closely with numerical simulations. Note that 
the discrepancy between values for Qrf in (a) and (c) is re- 
quired to fit the data, and likely arises from nonidealities in 
RF amplifier (Minicircuits ZHL-32A) frequency response in 
our experiments. For higher-order transitions 5 = nu, the 
missing resonances coincide with zeros of J§ / uj (2Qrf /^rf) 
as shown in the inset to Fig. 2a. The excellent agreement be- 
tween the frequencies of the observed missing resonances and 
theory confirms that they indeed arise from coherent destruc- 
tion of tunneling. 



FIG. 4: (a) A single trace from Fig. 2a at u RF = (27r)300kHz, 
shown with the expected quasistatic lineshape (red line) for Qr F = 
(2?r)15 MHz convolved with a Lorenzian of HWHM 0.5 MHz. (b) 
MW spectrum vs RF power. Solid lines illustrate functional depen- 
dence of Rabi frequency Qrf on RF power, (c) Resonant MW exci- 
tation (lumw ~ A) vs RF frequency, illustrating CDT. Experimental 
data obtained by averaging over the central 1 MHz line in data from 
Fig. 2a. Theory shows population in m s — after 5/xs numerical in- 
tegration of the Bloch equations for Eq. fusing Qmw — (27r)0.43 
MHz, Qrf = (27r)15.4 MHz incorporating 0.5 MHz pure dephas- 
ing. Upper curve indicates corresponding zeros in the Bessel func- 
tion J (2n RF /ujRF). 



driven with an effective Rabi frequency of ^mi^^ 




when the resonance condition S = uujrf is satisfied. The 
n-photon transition disappears for excitation frequencies ujrf 
given by the roots of the n^ order Bessel function, yielding 
the locations of the missing resonances. 

These missing resonances are related to so-called coherent 
destruction of tunneling (CDT), which is a strong-field phe- 
nomenon that only occurs when the Rabi frequency of an ex- 
citation exceeds its carrier frequency. CDT arises in systems 
obeying a Hamiltonian of the form of Eq. [2] when the off- 
diagonal transition (in this case Qmw) is small and the ampli- 
tude of the longitudinal oscillations (Qrf) is larger than the 
oscillation frequency (ujrf)- Physically, CDT can be inter- 
preted as a resonance in the Landau-Zener crossing between 
instantaneous eigenstates of Eq. [2] induced by a large ampli- 
tude oscillating detuning: it occurs when the amplitude and 
frequency of the oscillation conspire to precisely drive the sys- 
tem across the energy gap between instantaneous eigenstates 
during each half cycle. CDT was first theoretically recognized 
for a particle confined to a double well l28l ; it has since been 
observed in systems including cold atoms (l0l|29l and optical 
waveguides l30l . While recently there has been a theoretical 
description of CDT in NV centers 1311 . the phenomenon has 
not, to our knowledge, been previously observed in this con- 
text. 



III. CONSEQUENCES AND APPLICATIONS 

While observations of multifrequency excitation of the NV 
center present an opportunity to explore quantum optical phe- 
nomena, they may also have relevance to applications in quan- 
tum information science. In particular, our results highlight 
potential pitfalls in multiplexed electron-nuclear manipulation 
and point to possible applications using oriented microwave 
magnetic fields. Because the electron gyromagnetic ratio is 
several orders of magnitude larger than that of a nucleus, RF 
excitation intended for a nuclear spin will have a consider- 
ably larger effect on an electronic spin. Even with sequential 
pulses, care must be taken to control or account for phase ac- 
cumulation from the longitudinal components of the RF mag- 
netic field. Similarly, attempts to simultaneously manipulate 
electron and nuclear spins must account for CDT suppression 
of electron spin transition rates. 

Most spin resonance experiments rely on transverse excita- 
tion, but longitudinal coupling offers an additional technique 
for spin manipulation. In particular, it enables spin states to be 
brought quickly into and out of resonance without the need for 
fast magnetic field ramps. The longitudinal selection rules in- 
volve no spin flips, so this type of excitation could potentially 
be used to turn on and off a static transverse coupling. 

As an example relevant in the context of the NV center in 
diamond, the nonsecular terms in a contact hyperfine Hamil- 
tonian A(S+I- + S-I+) provide such a static transverse 
coupling between the electron-nuclear spin states |0,|) and 
| — 1, t). (Here, t, | encodes a nuclear spin 1 = 1/2 projec- 
tion along the NV axis, while numbers refer to the NV elec- 
tronic spin projection m s .) While a static mangetic field could 
be used to bring the initial and final states into resonance, 
leading to oscillations between them, a longitudinal magnetic 
field with strength Q cos ut oscillating at the energy differ- 
ence between these states will also lead to electron-nuclear 
spin flip-flop transitions at a rate ~ AQ/uo. With appropriate 
resonator ll32l or stripline lfT2l design, longitudinal excitation 
could thereby be used to implement electron-nuclear swaps or 
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extend optical nuclear spin polarization to arbitrary magnetic 
field strengths |7|. Although only one oscillating field is in- 
volved, this type of configuration mimics the multifrequency 
excitation discussed in this paper in that the nonsecular hyper- 
fine terms resemble detuned MW excitation (in the RWA) and 
the applied longitudinal excitation plays the role of the RF by 
enabling resonant transitions. 

IV. CONCLUSION 

Simultaneous application of weak MW and strong RF exci- 
tation in a two-level system leads to a rich set of phenomena 
beyond the usual rotating wave approximation. Using mul- 
tifrequency excitation of the NV center in diamond, we ob- 
serve quasi-static splitting, multiphoton resonances and co- 
herent destruction of tunneling in the electronic spin degree 
of freedom. Notably, the condition for strong-field behav- 
ior does not require Rabi frequencies comparable to the en- 
ergy separation between spin states, making this regime easy 



to access experimentally. Numerical simulations confirm our 
understanding that this behavior arises from combined lon- 
gitudinal and transverse excitation, with resulting selection 
rules allowing resonant two-photon transitions. These results 
indicate the importance of taking the orientation of applied 
microwaves into account when designing experiments; more- 
over, they invite exploration of control techniques based on 
longitudinal excitation. Ultimately, a greater understanding 
of multifrequency excitation may inform new techniques for 
manipulating spin systems. 
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